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Abstract 

The dynamics of solid deuterium (SD2) is studied by means of inelastic scattering (coherent and 
incoherent) of thermal and cold neutrons at different temperatures and para-ortho ratios. In this 
paper, the results for the generalized density of states (GDOS) are presented and discussed. The 
measurements were performed at the thermal neutron time-of-flight (TOF) instrument IN4 at ILL 
Grenoble and at the cold neutron TOF instrument TOFTOF at FRM II Garching. The GDOS 
comprises besides the hep phonon excitations of the SD2 the rotational transitions J = 1 
and J = 1 I— > 2. The intensities of these rotational excitations depend strongly on the ortho-D2 
molecule concentration Cq in SD2. Above E = 10 meV there are still strong excitations, which very 
likely may originate from higher energy damped optical phonons and multi-phonon contributions. 
A method for separating the one- and multi-phononon contributions to the density of states will 
be presented and discussed. 

PACS numbers: 28.20.Cz, 63.20.kk 
Keywords: neutron, phonons, solid deuterium 
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I. INTRODUCTION 



Solid deuterium (SD2) and solid hydrogen are typical quantum molecular solids. Each 
D2/H2 molecule exhibits large zero-point vibrations due to the small molecule mass. These 
quantum solids have been investigated theoretically il| and^by experimental techniques 
like inelastic neutron scattering j2| and Raman scattering ^]. Solid deuterium has an 
almost perfect hep crystal structure, when it is prepared under suitable conditions [4] (low 
pressure and T > 5 K). The phonon dispersion relation of solid ortho-deuterium (0-D2) was 
measured in the past by inelastic coherent neutron scattering [5|. These measured phonon 
dispersion relations are in good agreement with recent calculations based on molecular 

l]. In the case of solid para- hydrogen, the molecular 



dynamics ca. 
calculations 



iculations (MDC) 



l| disagree to some extent with the experimental data of Nielsen j^. The 
MDC predicts phonons above 10 meV in solid para-hydrogen, which did not appear in 
the neutron scattering measurements of Nielsen One possible explanation could be 
the restriction to energies below 9 meV in Nielsen's data analysis Recent inelastic 



incoherent neutron scattering data on solid para-hydrogen 
above 10 meV. 



^] are showing also excitations 



The D2 molecule has internal rotational modes, which are described by the rota- 
tional quantum number J. In the solid phase, J is still a good quantum number. Deuterium 
in states with even J ( J=0,2,4....) is termed ortho-deuterium (0-D2), whereas D2 in 

states with odd J (J=l,3,5, ) is termed para-deuterium (P-D2). At low temperatures 

(T ~ 6 K), about 99.999% of the deuterium molecules are in the ortho state, when thermal 
equilibrium is reached. At room temperature, D2 has an ortho concentration of 66.7%. 
After cooling down the D2 to the solid phase (T ~ 6 K), it normally takes months to reach 
the equilibrium of 99.999% 0-D2. This process can be accelerated by using paramagnetic 
materials like chromium oxides or hydrous ferric oxides 3, Isl, Iq] as catalysts for the 
para-ortho conversion. With this kind of converters, it is possible to reduce the conversion 
time to below one day in typical experiments and to adjust the 0-D2 concentration between 
66.7% and >= 98% for investigating in detail its influence on the structure and dynamics 
of solid deuterium. 
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Neutron scattering is an excellent tool to investigate the phonon system and the ro- 
tational transitions of D2 molecules in SD2. Neutrons are scattered by D coherently 



(^coh = 6.671 fm) and also incoherently {bi^c = 4.04 fm) [10[. Scattering of neutrons where 
a rotational transition (J \—>- J' , J ^ J') is involved, leads to an incoherent response of the 
system. The phonon excitations are present in the coherent part of the scattering, whereas 
the incoherent part is determined by rotational transitions and also by incoherent phonon 
scattering Q ■ 

In principle the cross section for neutron scattering in SD2 is described by, 

-j—bcoh Scoh (Q,^) + -j—binc Siac {Qi^)- (l) 



OUJOi I kg kg 

Here, ko is the wave number of the incoming neutron, k the wave vector of the scattered 
neutron, Si{Q,uj) {i =coh, inc) are the dynamical structure functions of coherent and 
incoherent scattering, respectively. The energy change of the neutron is E = Tiw , while 
hQ is the momentum transfer. The coherent and incoherent scattering lengths are 6coh and 
6inc respectiveley. 

In the incoherent approximation the neutron scattering cross section Eq. ([1]) can be 
expressed by 

The generalized density of states (GDOS) G{uj) comprises the complete phonon exci- 
tations of SD2 as well as rotational transitions of individual D2 molecules. Furthermore, 
multi-phonon excitations of the phonon system of SD2 jl4] should appear in the GDOS, as 
they are not corrected for in Eq. ([2]). The term e"^^*^'^-' is the well known Debye- Waller 
■actor, where W{Q) = \Q'^{u^), and (u^) ((m^) = o.25 A^) is the mean square displacement 
5I of the D2 molecule in the lattice. The quantity n{uj) is the Bose statistic function, and 
[n{uj) + 1] describes the creation of a boson, which causes the neutron energy loss. The 
scattering length &efr(Q) is a combination of coherent and incoherent scattering lengths. 
This combination depends strongly on the rotational transitions excited by the neutrons 
Q (see TABLE I). 
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Factorization (Eq. ([2])) is only applicable, if the total scattering response of solid 
deuterium contains an essential amount of incoherent scattering (incoherent approxima- 
tion), or if averaging over a sufficiently large Q range is performed. This is the case for 
thermal neutron scattering. 



TABLE I: Neutron scattering lengths b associated with rotational transitions [1/ 



J 


E [meV] 


62 


^ 1 


7.0 


0.375 62^, 


1^2 


13.5 


0.75 6L 


0^2 


21.0 


^coh + 0.6256L 


1^0 


7.0 


0.756L 



In principle, [bes{Q)] is expressed by 



(l-Cp 
+ (1 - Cp) ■ 



4 I + -bt 



foiQas/2) 



-bL{jl{Qas/2)) 



(3) 



+ Cp ■ (2j?(gas/2) + 3jl{Qas/2))] 

+ c, ■ + €c) (jo(Q«./2) + 2j1(Qa,/2))] , 

where ji{Q ■ as/2) are the spherical Bessel functions of order i. The first term in 
Eq. (IH]) describes the transition J = (even-even), while the second term the 
J = (— i> 1 (even-odd) transition. The third and fourth term take the J = 1 2 (odd-even) 
and J = 1 1-^ 1 (odd-odd) transition into account. The parameter = 0.74 A is the 
distance of the deuterons within the D2 molecule [5], while Cp is the concentration of the 
molecules in the para state {cp + Co = 1, Co is the concentration of molecules in the ortho 
state). 



With the measured inelastic neutron cross section and Eq. ([2]), it is possible to de- 
termine the GDOS. The GDOS for differently prepared SD2 crystals with different ortho 
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concentrations are the main subject of this paper. Furthermore, a methode of extracting 
the density of states of one-quasi-particle excitations Gi{E) from our data is presented. A 
detailed analysis of the dynamical and static structure of solid deuterium will be presented 
in a forthcoming publication. 



II. EXPERIMENTAL DETAILS 



The experiments on inelastic neutron scattering were performed at the Time-Of-Flight 
(TOP) spectrometer for thermal neutrons IN4 15|] of the Institute Laue-Langevin (ILL), 
Grenoble, France and at the cold TOP spectrometer TOFTOF [L^ at FRM II, Garching, 
Germany. 



The measurements at IN4 were carried out at two different wavelengths of the in- 
coming neutrons (A ~ 2.2 A and 1.1 A). The energy resolution for the two different 
wavelengths had been determined to be AE2.2 = 0.7 meV and AEi ^ = 3.4 meV, using the 
standard technique of elastic neutron scattering on a vanadium sample. 

The IN4 spectrometer supplies a typical thermal neutron flux of 5 ■ 10^ cm^^ s^^ on 
the sample. The beam size on the sample is ~ 2x4 cm^. The scattering angle 26 varies 
between 13° and 120°. The beam divergence is A^^ = 1°. 



The sample cell for the solid deuterium was designed and constructed to fit into the 
standard "orange cryostat" of ILL [l7|. The sample cell is cooled down to the desired 
temperature by using an exchange gas (helium) in the chamber which surrounds the 
sample cell. The cooling inside the cryostat is done by a liquid-nitrogen shield and a 
liquid-helium heat exchanger. Neutron scattering on solid/liquid deuterium was performed 
at temperatures between T = 3.5 K and T = 23 K. The temperature stability was better 
than 0.1 K. 

The SD2 is frozen inside the double- wall cylinder volume (see Fig. [1]). The thickness 
of the SD2 in the cell is fixed to 5 = 2 mm. The attenuation of the incoming thermal 
neutrons {Eq ~ 17 meV) in the filled sample cell is ~ 25%, and therefore multiple neutron 
scattering inside the sample can not be neglected. The influence of multiple scattering will 
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be discussed in detail in the section " Results and discussion" . The sample cell is connected 
to the D2 gas handling system via a 12 mmxl mm stainless steel pipe. 



The gas handling system is especially designed to keep the pressure of the D2 gas flow into 
the sample cell quasi-independent of the temperature within the cell. This is done by using 
needle valves for the gas flow. With this technique, it is possible to follow a special path in 
the phase diagram (see Fig. [2]) of D2. This can be either the path from the gas via the liquid 
to the solid phase, or the path directly from the gas to the solid phase (sublimation). The 
triple point of D2 is at [Tt = 18.7 K, Pt ~ 180 mbar]. At pressures below Pt, direct sublima- 
tion of D2 is possible, while at pressures above Pt, the D2 is liquefled before it becomes a solid. 

The para-to-ortho converter j^, which accelerates the conversion from the P-D2 to 
0-D2 state (Fig. [3]), is attached to the gas handling system via a stainless-steel pipe. This 
converter is an independent cryogenic system, cooled by a cold flnger. The main part of 
this system is a copper cup, fllled with a paramagnetic powder. This powder was extracted 
from an OXISORB [9] unit, which is normally a gas cleaning tool. OXISORB contains 
paramagnetic chromium composition, which acts inside the converter as a catalyst. The 
upper part of the copper cup is closed by a sinter disc, which keeps the powder in the 
cup, but is permeable to D2 gas. This copper cup is mounted on the flnger of the cooling 
machine, and can be operated between T = 10 K and room temperature. During the flUing 
mode, the cup is fed with D2 gas from the gas handling system. The D2 is liquifled inside 
the cup, and is kept there at the boiling point for the necessary conversion time (hours), 
keeping the liquid at the boiling point shortens the conversion time considerably (factor 2) 



The fllling procedure of the sample cell was monitored with the help of the IN4 
spectrometer. Filling was stopped when the scattered-neutron intensity reached its 
maximum and stayed constant. This method was used for fllling natural deuterium gas 
from the dump and also for fllling the converted deuterium gas from the converter unit. 

The experiments at TOFTOF of FRM II were performed with an equipment similar 
to that used at IN4. These experiments were performed before the IN4 experiment and 
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served to optimize the sample cell for the IN4 experiment. Sample cell and cooling machine 
at the TOFTOF experiment were different from the setup at IN4. The thickness of the 
SD2 in the cell was bigger than that (5 = 3 mm) in the IN4 sample cell. The attenuation 
of the incoming thermal neutrons {Eq ~ 20 meV, A ~ 2.0 A) in the filled sample cell is 
~ 30%, and therefore multiple neutron scattering inside the sample can not be neglected. 
The influence of multiple scattering will be discussed in detail in the section "Results and 
discussion". Later, after the IN4 experiment, a second experiment was done with SD2 at 
the TOFTOF, in order to improve the statistics of neutron data for energy transfers larger 
than 10 meV. 

The sample cell at TOFTOF is cooled by a closed cycle cold-head machine with He exchange 
gas. With this setup temperatures down to T ~7 K were reached. The measurements 
at TOFTOF were carried out at three different wavelengths of the incoming neutrons 
(A ~ 2.0 A , A ^ 2.6 A and A ^ 6.0 A). The energy resolution for the different wavelengths 
has been determined to AE2.0 = 1.23 meV, AE2.G = 0.74 meV and AEq,q = 0.07 meV, 
using the standard technique of elastic neutron scattering on a vanadium sample. 
At TOFTOF the typical cold neutron flux on the sample amounts to 7 ■ 10^ cm~^ s^^. The 
beam size on the sample is ~ 3 x 5 cm^. The TOFTOF detector bank covers a range of 
26* = 7.5° ... 140°. The beam divergency is A6' ~ 1°. 



III. SAMPLE PREPARATION 

The structure and dynamics of SD2 may be influenced by the way of preparing the 
crystals fisl. [lol. Therefore samples of different character were prepared by different freezing 
procedures. Our own experiments had shown differences in the optical quality of D2 for 
different freezing procedures (see Fig. H]). 

Freezing experiments with p-sH2 and 0-SD2 films on MgF2 coated sapphire plates jiol 
showed very interesting features concerning the microscopic and mesoscopic structure of 
solid deuterium and hydrogen films. Below < 6 K (deposition temperature) a mixture 
of hep and fee structures is possible in the SD2 film. The fee component disappears 
if the temperature is raised above 9.5 K (annealing). The remaining hep structure is 
preserved, even if the temperature is lowered again below 9.5 K. These observations were 
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one motivation to look for a possible influence of the freezing procedure on the structure 
and dynamics of SD2. The following freezing procedures were applied: 



Liquid-Solid (LS): D2 is liquefied in the sample cell at T ~21 K. The filling is mon- 
itored by neutron scattering measurements at IN4. After filling, the temperature of the 
cooling controller is lowered to a temperature slightly below T = 18.7 K and kept constant, 
the sample starts to freeze slowly. The SD2 sample is then cooled down slowly to the desired 
measurement temperature. Slow cooling of the SD2 is necessary to avoid sudden cracking of 
the solid deuterium. 

Liquid-Solid-Melt-Solid (LSMS) : An alternative method of producing good samples 
of SD2 is remelting the solid at T = 18.7 K. When the solid is melting, a mixture of liquid 
and small pieces of solid deuterium crystals will emerge. The character of this mixture can 
be monitored by neutron-scattering measurements. The static structure factor indicates the 
crossover from sohd/liquid to complete liquid. The basic idea of this method is to stop the 
melting shortly before all solid pieces of SD2 disappear. The remaining solid fragments serve 
as crystallization seeds for a new solid-deuterium sample. 

Gas-Solid (GS): The re-sublimation of deuterium from gas is done be keeping the 
sample-cell temperature at T ~ 7 — 9 K and maintaining a constant gas flow. With this 
method it is possible to grow slowly a solid-deuterium sample in the cell. How far the cell 
is filled with solid material is again monitored by neutron-scattering measurements. 

Turbo Solid (TS): Fast freezing of SD2 from the liquid phase at a temperature below 
T ~ 10 K is leading to an opaque solid-deuterium sample (see Fig. Hb). The analysis of the 
static structure, extracted from the neutron scattering data, shows an almost ideal powder 
diffraction pattern for this "turbo solid". A detailed analysis of the structure of SD2 will be 
presented in a forthcoming paper. 
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IV. RESULTS AND DISCUSSION 



The TOF data were corrected for the empty-cell measurements, and normalized to the 
incoming flux and the vanadium standard. Based on this data the GDOS were determined 
by the IN4 computer code gdos. This code is implemented in the data analysis package 
LAMP 20|, which is used for treating the data obtained from neutron scattering experiments 
at ILL. The basic method of this calculation is a sampling over a large Q range (neutron 
energy loss side), using the incoherent approximation |21j] and referring to Eq. ([2]). 
The GDOS is calculated by: 

SiQi< 26 >),u;) 
GDOS{.) = . . (4) 



The term Q{< 26 >) is a Q value for each energy transfer hw considering an averaged 
scattering angle < 26 >. This simple approximation for the GDOS is commonly not equal 
to the real density of states (DOS). In many cases the agreement between the GDOS and 
DOS is remarkable. The GDOS incorporates the one-phonon and multi-phonon excitations. 
In the flrst order it is not so easy to separate the one-phonon part from the multi-phonon 
part. Nevertheless the GDOS is a useful tool to investigate the spectrum of excitations of 
the sample. 

The experimental values of the GDOS for natural deuterium (0-D2 66.7%) in the liquid 
and also in the solid phase for different preparation methods of the crystals are shown in 
Fig. [5l The data for almost pure (co = 95%) 0-D2 are shown in Fig. [6l 

The GDOS for differently prepared crystals at the same 0-D2 concentration does not 
show variations, beyond statistical effects. It seems that the way of preparing the crystals 
does not change the spectra of phonon excitations in SD2. The main difference between 
the GDOS for different 0-D2 concentrations is the strength of the rotational transition 
J = I— s> 1. Increasing the 0-D2 concentration enhances this transition. The GDOS shows 
for both 0-D2 concentrations a peak at E =12. .13 meV. Furthermore, at Co = 66.7% there 
is a hint for the rotational transition J = 1 2. This excitation is not visible for Co = 95%. 



An example for the influence of the 0-D2 concentration Cq on the rotational transitions is 
shown in Fig. [71 These data are taken from the later TOFTOF measurements (A ~ 2.0 A). 
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It is evident that increasing Co leads to a larger cross section for the J = (—> 1 transition 
and also to a decrease of the J = 1 i-^ transition. 

The main focus of the data analysis presented here is the extraction of the density of states 
(DOS) Gi{E) of phonons and rotational transitions in solid D2 ("quasi particle picture"). 
The basic idea is to compare the integral over the scattering angle of the neutron scattering 
data da/dEf with a calculated neutron cross section, using the incoherent approximation 



developed by Turchin [2^ for a cubic crystal with one atom in the primitive cell. With the 
aid of this theory, it is possible to calculate the one-phonon and multi-phonon contributions 
to the neutron cross section da/dEf, if Gi{E) is known. 

This approach is only valid, if the detected neutrons are scattered neutrons from all possible 
directions of the first Brillouin zone in the crystal. This is normally the case, if one has 
a powder sample of the crystals. The diffraction pattern of fast frozen deuterium (TS, 
Cp=33.3%, T=9.5K) shows an powder like behavior, although some small texture effects are 



visible (see Fig. [8]). A Rietveld fit on this diffraction pattern gives the right values 2^ for 
a and c of a hep solid deuterium crystal (a =3.596 +/- 0.005 A, c =5.860 +/- 0.009 A). 
The positions of the different Bragg peaks correspond to a hep solid D2 crystal. The fit- 
model contains also the Debye- Waller factor, because at larger scattering angles (higher 
Q-values) the scattering intensity is significantly reduced. The fit result for the mean square 
displacement < >=0.208 +/- 0.011 is smaller compared to value reported by Nielsen 
j^l. These results induced the decision to use only data from TS samples (second TOFTOF 
measurement) to determine the density of states Gi{E). 



00 

as{Ei ^Ef)=Y, (yn{E, Ef), (5) 

n=l 



ao(Q)=47r.[6eff(Q)]' (7) 

The summation over the index n (number of quasi particles involved in the scattering) 
contains all higher-order multi-quasi-particle excitations like multiphonons. With Ei and Ef 
are the energy of the incoming and outgoing neutron in the scattering process, respectively. 
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fi is the mass number of the D2 molecule (/i = 4) and e is the energy transfer {E = Ei — Ef) 
to the scattered neutron. 



/+00 
U_,{E')f{E-E')dE', (8) 
-00 

f(E) = ^^^'^'^ (9) 

The function f{E) (see (Eq. ([9])) contains the density of states Gi{E), and it is possible 
to calculate the higher order functions fn{E) with the help of Eq. ([8]). 



Fn{Ei,Ej)= ^ / x"e-^/"da; (10) 



In Eq. ( fTOl) Fn{Ei, Ef) is a kinematic factor, resulting from the integration over all angles 
in the neutron scattering process. The quantity r is the average energy of the excitations 
(mainly phonons) in SD2. 

The starting point of the determination of Gi{E) is a first-guess model for Gi{E), which 
should contain the phonons as well as the rotational transition J \—>- J'. Following the 
approach of Yu et al. [2^ a model with seven Gaussian functions is used to parameterize 
Gi{E). The contribution of the three-particle excitations {da^/dEf) was approximated by 
an additional Gaussian function. 

Gi(i?) = J]«,-=^e (12) 

Treating phonons and rotational transitions as independent excitations may be not 
strictly correct in the calculation of the multiphonon contribution, but an eventual cor- 
rection would occur at the far end of the energy spectrum above E > 14 meV, which is not 
crucial important here. 

The neutron cross section da/dEf is calculated with Eqs. ([5]) - f|T2l) and compared with 
the measured NQ-dau/dEf {Nq is a normalization factor). The parameters [ai,ai,ei;l = 1..7] 
are determined by minimizing the difference between measured and calculated cross section. 



12 



Multiple scattering is in the case of the TOFTOF measurement not negligible. The deu- 
terium sample scatters approximately 25% of the incoming beam. The scattering intensity 
is 50% elastic and 50% inelastic. The contribution of multiple scattering is a fiat falling 
background {E >3 meV) and peak-like at E ~0 meV. The multiple scattering was calcu- 
lated (convolution of the scattering model with itself) by using the experimental deduced 
Gi{E) and applying the Turchin theory for one- and two-phonon neutron cross section. 
The magnitude ( 13%) of the multiple scattering was determined by using the work of V.F. 



Sears 



25l | ("rule of thumb" 26|]). The flat falling background is smooth and has almost no 



texture. This background was included in the fit model, and cross checked afterwards with 
the convolution of da/dE with itself. The magnitude of multiple scattering, deduced from 
the fit must be in the same order, as predicted by the "rule of thumb" -estimation. 

The TOFTOF measurements (see Fig. [9]) clearly show the onset of two-particle contribu- 
tions above ~ 5 meV and also one-particle excitations above i^^ = 10 meV. The TOFTOF 
data indicates three-particle excitations above E = 12 meV. The peak at -E = 5 meV ob- 
viously originates from TA-phonons [2,], whereas the peak at ~ 7.5 meV is a signature of 
the rotational transition J = i— * 1. 

Figure [To] shows the result for the one-particle density of states for solid D2 at Co = 66.7% 
and Co = 98% (area normalized to 1). Part a) of Fig. [10] shows the convoluted Gi{E), 
while in part b) Gi{E) is de-convoluted with an Gaussian pseudo-resolution function with 
an FWHM =1.0 meV. This value was deduced from the FWHM of the elastic peak of 
da/dE. 

The influence of the ortho concentration Co on the density of states Gi{E) manifests itself 
(Fig. [To]) in the expected enhancement of the peak at ~ 7.5 meV for larger Cq. Increasing 
the number of ortho molecules in solid D2 leads to a higher probability for J = ^-^ 1 
transition. The position of the J = 1— * 1 transition depends slightly on the concentration 
of para molecules Cp in the SD2. A larger Cp lowers the J = ^ 1 energy {6E ~ 0.3 meV for 
Cp = 33.3%). This energy shift can be explained by the interaction between P-D2 (J = 1) 
and 0-D2 molecules (J = 0) [4]. 

The DOS derived from our measurements differs remarkably from earlier published results 



2J], which used the phonon dispersion measurements of Nielsen [5] to calculate the DOS 



of solid D2. The DOS published by Yu et al. [24] contains pronounced peaks at 5 meV and 
9 meV, which were identified with the acoustic and optical phonons in the hep crystal. The 
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J = 1-^ 1 transition was not considered by Yu et al. [2^, and the DOS vanishes above 
10 meV. The optical phonon contribution of the actually measured DOS for Cp=33.3% 
seems to be smeared-out in the region of 8 — 12 meV, and does not peak at ~ 9 meV, as 
it was published by Yu et al. [24,]. This behavior was already observed in solid hydrogen 
by incoherent neutron scattering of Bickermann et al. 27| and Colognesi et al. 



Bickermann et al. explain the smearing-out of the higher-energy optical phonon groups 
by the anharmonicity of the quantum- crystal solid H2. This anharmonicity leads to larger 
phonon line widths at higher energies. The DOS of Colognesi et al. shows at ~ 9 — 10 meV 
a small bumb, which could be identified as optic phonons in solid hydrogen. In the case of 
almost pure 0-D2 (Cp=2%) a peak at ~ 9 meV in the DOS was found in the fit of our 
data. The amplitude of this optical peak is smaller than as it was published by Yu et al. 
The appearance of multiphonons above 5 meV in our measurements leads to considerable 
line widths of the phonon groups above this certain energy. Nielsen's DOS was determined 
by a pseudoharmonic theory, which takes care of the renormalization of the phonon energies, 
but does not consider the change of the phonon line widths 27| |. 

A second new feature of the result for Gi{E) presented here is a clear peak-like structure 
at ~ 12 meV. The intensity of this peak is more pronounced at higher Cq. This peak 
might be caused by the excitation of a combination of the rotational transition J = 1— > 1 
(~ 7.5 meV) and phonons with ~ 5 meV energy. The DOS for phonons with ~ 5 meV is 
very large. A combined excitation of phonons and rotational transitions of the D2 molecules 
is certainly a result of the interaction of the rotational transitions with the lattice of the 
solid D2 crystal (phonon-rotation coupling) 29!]. A higher energy phonon {E > 12 meV) is 
able to excite a low energy phonon and also the rotational transition J = ^— 1 through 
the phonon-rotation coupling. This decay of higher-energy phonons was maybe seen in our 
data. This peak could be on the other hand maybe an artefact, resulting from the incomplete 
separation of the one-particle-excitations and two-particle-excitations in the DOS, and has 
to be further investigated. 

The DOS of sohd D2 with higher para concentration (cp =~ 33.3%, Co ~ 66.7%) contains a 
peak structure (see Fig. [TOl) at ~ 14 meV, which is originated by the rotational transition 
J = 1 ^— s> 2. This peak is not seen at Cp = 2%, Cq = 98%. This effect is obvious, because the 
J = 1 2 transition scales with the number of para molecules in the solid D2. 
The result for Gi{E) derived from our neutron scattering data can be used to calculate 
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the mean square displacement < >. This value (< >=0.22 for Co ~ 66.7%) is 
quite close to the result from the Rietveld fit on the neutron diffraction data of TS samples 

n 

(co ~ 66.7%, < >=0.208 A^), and smaller than the result reported by Nielsen 5|]. In 
this context it is interesting to mention the recent published work of Bafile et al. 30||, who 
reports on neutron diffraction measurements of solid deuterium close to melting. The mean 
square displacement increases to < >=0.43 at T ~18.7 K. The conclusion of Bafile 
et al. is the break down of the pseudoharmonic approach for solid deuterium close to the 
melting point. 



V. CONCLUSION 



In this manuscript, we reported investigations on the density of states (DOS) of solid 
D2 for different 0-D2 concentrations with the powerful tool of neutron scattering. The 
experiments were performed at time-of- flight instruments at ILL, Grenoble and at FRM 
II, Garching. The DOS was extracted from the neutron data by applying the incoherent 
approximation. The results for the GDOS (generalized density of states) and the DOS 
show interesting details, which have not been published previously. One feature is the 
appearance of the J = 1— > 1 in the DOS. The strength of this excitation increases with 
the 0-D2 concentration Cq. The DOS of solid D2 contains a strong phonon signal at 5 meV, 
but weak signals of optical phonons. Our results for the DOS indicate a smearing-out of 
the higher energy optical phonons for higher P-D2 concentrations, which was also reported 
for solid hydrogen with high 0-H2 concentrations. The reason for this effect may be 
explained by the anharmonicity of the solid D2 crystal. Furthermore we see two-phonon 
processes at considerably low energies in our data (onset at 5 meV). Another new result 
is the appearance of a peak at ~ 12 meV in the DOS of solid D2. This peak increases 
with higher 0-D2 concentration, and could be explained by a coupled excitation of a 
phonon and a rotational transition, or by an possible incomplete separation of the one- and 
two-particle-contribution. This effect has to be further investigated. Our data of natural 
deuterium (co ~ 66.7%) are showing the rotational transition J = 1 ^ 2. 
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Figure captions: 

Fig.[U 

Sketch of the D2 sample celL 



Fig.H 

D2 phase diagram with different process paths studied: a) hquefaction from the gas, b) 
sohdification from hquid, c) subhmation at high temperature, d) direct condensation. 



Fig. El 

Sketch of the para-to-ortho converter [9|]. The copper converter cup is filled with OXISORB 
grain, and mounted on the cold-finger of the cooling machine. 



Fig.H 

Pictures of solid deuterium samples for different freezing procedures. 

Crystals from the liquid phase are transparent, if they are frozen slowly, while a fast 
freezing/cool down (within minutes) leads to a non-transparent solid. 

a) solid deuterium slowly frozen from the liquid phase (LS), 

b) solid deuterium slowly frozen from the gas phase (GS), 

c) solid deuterium fast freezing/cool down from the liquid phase (TS). 



Fig.EJ 

Generalized density of states (GDOS) of natural deuterium (cq = 66.7%): 

liquid D2 (□) at T =21 K, solid D2 (A), rapidly frozen out from the liquid phase and fast 

cooled down, solid D2 (O)) slowly frozen out from the liquid phase, solid D2 (^), slowly frozen 

out from the gas phase. All solids have a temperature of T =4 K. Peaks are convoluted with 

IN4 energy resolution curve: Normalized by GDOS{E) ■ dE = 1. 

Each GDOS is separated by a shift of 0.1. 

Fig. El 

Generalized density of states (GDOS) of converted deuterium (co = 95%): 
liquid D2 (□) at T =20 K, solid D2 (A), slowly frozen out from the liquid phase, solid D2 
(O)) slowly frozen out from the gas phase, solid D2 (*), frozen out at the melting point and 
slowly cooled down. All solids have a temperature of T =4 K. Peaks are convoluted with 
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IN4 energy resolution curve: Normalized by GDOS{E) ■ dE = 1. 

Each GDOS is separated by a shift of 0.1. 

Fig.B 

Dynamical neutron scattering cross section of solid D2 for Cq = 66.7% (□) and Cq = 98% 
(O) at T = 7 K. Data from the TOFTOF measurements. 



Fig. El 

Neutron diffraction pattern data of fast frozen solid D2 (TS) for Cq = 66.7% (A) at 
T = 9.5 K and comparison with a Rietveld fit (solid line) for powder like hep solid D2. 
Data from the TOFTOF measurements. 



Fig.l 

Dynamic neutron scattering cross section of solid D2 for Cq = 66.7% and Cq = 98% at 
T = 7 K. 

Comparison of data with calculated neutron cross sections. 

The one-particle contribution is shown by the dashed line, the two-particle contribution by 
the dotted line, and the three-particle contribution by the dash-dotted line. Contribution 
of multiple scattering is shown by the dot-dot-dash line. Data from the TOFTOF 
measurements. 



Fig.[ini 

a) One-particle density of states of solid D2 for Cq = 66.7% (O); Co = 98% (□) at T = 7 K. 

b) Comparison of de-convoluted DOS with data {-k) from Yu et al. 24 1. 

The DOS (part b)) are de-convoluted with the FWHM of the elastic peak of da/dE. 
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FIG. 1: 

Sketch of the D2 sample cell. 
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FIG. 2: 

D2 phase diagram with different process paths studied: a) hquefaction from the gas, b) 
sohdification from hquid, c) subhmation at high temperature, d) direct condensation. 
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Sketch of the para-to-ortho converter 



Q. 



FIG. 3: 



The copper converter cup is filled with OXISORB 



grain, and mounted on the cold-finger of the cooling machine. 
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FIG. 4: 

Pictures of solid deuterium samples for different freezing procedures. 
Crystals from the liquid phase are transparent, if they are frozen slowly, while a fast freezing/cool 
down (within minutes) leads to a non-transparent solid, 
a) solid deuterium slowly frozen from the liquid phase (LS), 
b) solid deuterium slowly frozen from the gas phase (GS), 
c) solid deuterium fast freezing/cool down from the liquid phase (TS). 
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FIG. 5: 

Generalized density of states (GDOS) of natural deuterium (co = 66.7%): 
liquid D2 (□) at T =21 K, solid D2 (A), rapidly frozen out from the liquid phase and fast cooled 
down, solid D2 (O); slowly frozen out from the liquid phase, solid D2 slowly frozen out from 
the gas phase. All solids have a temperature of T =4 K. Peaks are convoluted with the IN4 
energy resolution curve: Normalized by GDOS{E) ■ dE = 1. 
Each GDOS is separated by a shift of 0.1. 
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FIG. 6: 

Generalized density of states (GDOS) of converted deuterium (co = 95%): 
liquid D2 (□) at T =20 K, solid D2 (A), slowly frozen out from the liquid phase, solid D2 (O)) 
slowly frozen out from the gas phase, solid D2 (*), frozen out at the melting point and slowly 
cooled down. All solids have a temperature of T =4 K. Peaks are convoluted with IN4 energy 
resolution curve: Normalized by GDOS{E) ■ dE = 1. 
Each GDOS is separated by a shift of 0.1. 
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FIG. 7: 

Dynamical neutron scattering cross section of solid D2 for Co = 66.7% (□) and Co = 98% (Q) at 

T = 7 K. Data from the TOFTOF measurements. 
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FIG. 8: 

Neutron diffraction pattern data of fast frozen solid D2 (TS) for Cq = 66.7% (A) at T = 9.5 K 
and comparison with a Rietveld fit (solid line) for powder like hep solid D2. Data from the 

TOFTOF measurements. 
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FIG. 9: 

Dynamic neutron scattering cross section of solid D2 for Co = 66.7% and Co = 98% at T = 7 K. 

Comparison of data with calculated neutron cross sections. 
The one-particle contribution is shown by the dashed line, the two-particle contribution by the 
dotted line, and the three-particle contribution by the dash-dotted line. Contribution of multiple 
scattering is shown by the dot-dot-dash line. Data from the TOFTOF measurements. 
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FIG. 10: 

a) One-particle density of states of solid D2 for Cq = 66.7% (O); Co = 98% (□) at T = 7 K. 
b) Comparison of de-convoluted DOS with data (*) from Yu et al. [24 1. 



The DOS (part b)) are de-convoluted with the FWHM of the elastic peak of da/dE. 
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